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Abstract. We propose a method for modeling the magnetic properties of nanomaterials with 
different structures. The method is based on the Ising model and the approximation of the 
random field interaction. It is shown that in this approximation, the magnetization of the 
nanocrystal depends only on the number of nearest neighbors of the lattice atoms and the values 
of exchange integrals between them. This gives a good algorithmic problem of calculating the 
magnetization of any nano-object, whether it is ultrathin film or nanoparticle of any shape and 
structure, managing only a rule of selection of nearest neighbors. By setting different values of 
exchange integrals, it is easy to describe ferromagnets, antiferromagnets, and ferrimagnets in a 
unified formalism. Having obtained the magnetization curve of the sample it is possible to find 
the Curie temperature as a function of, for example, the thickness of ultrathin film. Afterwards 
one can obtain the numerical values for critical exponents of the phase transition "ferromagnet 
- paramagnet" . Good agreement between the results of calculations and the experimental data 
proves the correctness of the method. 



1. Introduction 

Creation of magnetic materials with predetermined properties - the task of unquestionable 
importance. While the manipulation of individual atoms becomes a routine, nanotechnologies 
need appropriate calculation methods, which would allow prototyping magnetic properties of 
materials without heavy calculations on supercomputers. The method proposed below satisfies 
such requirements. Moreover, it allows calculus of the properties of magnets of any type 
(ferromagnets, antiferromagnets, ferrimagnets) to be provided using a uniform formalism. 

Another feature of nano-objects is the existence of size effects at their scale. Obviously, such 
effects are different in two-dimensional films and three-dimensional particles. Furthermore, it 
is common to use different methods of modeling for nanoparticles of different shape. But if to 
consider the theory with short-range interaction, the global geometry of the material should not 
play any role. In this case, there must be a simulation method that is insensitive to the shape 
of the simulated sample as a whole. We offer in this paper such a method. 

2. General model of magnetic material 

We assume that magnetic atoms are distributed over N sites of the sample with probability p. 
According to [T], the distribution function for random interaction fields H on a particle located 



at the origin can be defined as: 



W (#) = Js(^H-Y^h k (r k , m fc ) fj ^ « ( r fc " r*o)J * fe dm A , (1) 

where 5(x — Xq) - Dirac delta function, h k = % (m^, r^) - field created by atoms with 
magnetic moments located at coordinates r^, r k o _ the coordinates of lattice sites, F(m k ) - 
the distribution function for the magnetic moments, which in the approximation of Ising model 
for a ferromagnets can be represented as follows: 

F (m fc ) = (a k 5 (6 k ) + fa8 (9 k - vr)) ((1 - p) 5 (m k ) + p6 (m k - mo)). (2) 

Here k - the angle between m k and OZ-axis, a k and f3 k - relative probabilities of the spin 
orientation along and against OZ-axis (9 k = and 9 k = it, respectively); mo - magnitude of 
magnetic moment of a magnetic atom. Probabilities a k and f3 k hold normalization condition 
ot k + Pk = 1- I n the approximation of nearest neighbors and the direct exchange interaction 
between magnetic atoms, the equation can be represented as: 
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W j (H) = Y,P z ~ n (l-p) n E E "iJiH-M^Jj), (3) 

where kj is a set of nearest neighbors of the magnetic atom numbered as j, z = dimfcj - its 
coordination number; C" (kj) - subset of n atoms of the total number of z nearest neighbors 
of j th atom; L (Q) is a binomial set of permutations of an arbitrary set £1 with the amount 
of elements equal to 2 dimf1 . Introducing symmetric notation a_ n = 1 — a n we have got 
= Y[ V {v=v,l v £ki) a ±lv and M h = T,nek 3 ±m n = rnoEn&k, .±|2a n - 1|. Finally, Jj is the 
constant of exchange interaction (exchange integral). 

Using the expression for the distribution function of interaction fields ([3]), one can obtain 
equations that determine average relative magnetic moments at each site of lattice: 



ft^(^)w j (H)dH = j2p z - n a-p) n E E ^ tanh (T^} 
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Expression Q allows to investigate the dependence of total magnetic moment of the sample 
M = Y2j=i a ^ t ne temperature T and concentration p, as well as to determine the dependence 
on the number N of atoms of the temperature of phase transition and the percolation threshold. 
System of N equations with N unknowns Q can be solved numerically using Newton's methocQ 

3. Modeling of magnetic materials with different crystalline structure 

Equation Q, despite the complicated form, has written in the algorithmically convenient form. 
This form allows to simulate the magnetic properties of materials, based only on the knowledge 
of the crystalline structure of the sample and the numerical values of the exchange integrals. 

Crucial part of this method is the rule of selection of nearest neighbors for each atom (the 
way of constructing of the set kj). By changing this rule of selection we can easily adjust our 
model to different physical systems. We also can simulate antiferromagnetic and ferrimagnetic 
materials by reversing signs and values of exchange integrals of individual atoms. 

Consider two examples of applying the method described above. 

1 In our research we use Wolfram Mathematica for fast lattice prototyping and C++ and Fortran non-linear 
solvers with python wrapper for precise calculus. 
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3.1. Ultrafine particles 

Consider a nanoparticle of 7Y atoms. Then Q is a system of ./V independent equations with 
7Y unknowns. In some cases, when p = 1, the symmetry of the particle reduces the number of 
unknowns, but in general all variables are different. 

For example, for cubic-shaped nanoprticle with n atoms on the edge and simple cubic lattice 
there is a system of n 3 non-algebraic equations with n 3 unknowns. Even relatively small cubic 
particle with 10 atoms on the edge contains 10 3 atoms, which makes calculation non-trivial and 
significantly reduces its accuracy. 

3.2. Thin films 

Now consider ultrathin film that is composed of N infinite monolayers. In a simple case when all 
atoms from the same layer are equal, Q turns into a system of N equations with N unknowns. 
Study of size effects in films is much easier than in the particles, because it is sufficient to 
simulate a film of 10-15 layers to achieve the bulk properties. (See fig. [TJ) 

General equation that determines the average relative magnetic moment fi n in n-th monolayer 
is given by Q. Replacing in expressions for uj£ and Mf all a± n on their average values 
(a± n ) = (1 ± /U n )/2 and substituting ^ into Q, one can obtain the equations that determine 
fj, n in each monolayer: 
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where z n ,n is the number of nearest neighbors in the n-th layer, z n -\^ n is the number of nearest 
neighbors of the atom in (n — l)-th layer, located in the n-th layer; i nn = JnnWn/ 1 J\\m\-, 
in—l,n = Jn— i,nm n -i/Jiimi, 

i"n,n+l — Jn, n+im n +i/ Jumi, t = kT/Jixmi. Using ([5]), one 
can study the dependence of the average magnetic moment of the film on its temperature and 
thickness [2], and, as a consequence, the dependence of Curie temperature on the thickness (see 
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Figure 2. The dependence of reduced 
Curie temperature Tc on thickness 
of ultrathin film (JV is the number 
of monolayers) for different crystal 
lattices (FCC, BCC and SC) and 
different crystallographic orientations 
of the surface. 



4. Comparison with experimental data 

Comparison of the results of modeling with experimental data was carried out directly and 
indirectly. 

4-1- Direct comparison: relative Curie temperature change 

In P], the dependence of Curie temperature Tc(N) on the size N of cubic-shaped particle 
was calculated. And then the magnitude of the relative change of the Curie temperature 
e(N) = Tc(N)/Tc(N — > oo) was obtained. Comparison with experimental data [1| revealed the 
validity of the modeling method. 

4-2. Indirect comparison: critical exponents 

It is possible to obtain the value of critical exponent v of spin-spin correlation for phase transition 
"ferromagnet - paramagnet" by appropriate approximation of e{N) defined above. Critical 
exponents for ultrathin films of different crystalline structures were obtained in [2]. It was 
shown that the value of v for three-dimensional Ising model is independent from the type of the 
lattice. The numerical value of v is close to the value obtained from RG-calculations [5]. 

5. Conclusion 

Offered method of simulating allows to substitute the real magnetic material with the relatively 
simple model, where the key role is played by the coordination number and the exchange 
integrals. All the magnetic geometry lies in these two integral characteristics. By setting only 
the rules of choice of nearest neighbors and values of corresponding exchange integrals, we make 
the modeling procedure very effective. 
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